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For the baseline design of the advanced Laser Interferometer Gravitational-wave Observatory 
(LIGO), use of optical cavities with non-spherical mirrors supporting flat-top ("mesa") beams, po- 
tentially capable of mitigating the thermal noise of the mirrors, has recently drawn a considerable 
attention. To reduce the severe tilt-instability problems affecting the originally conceived nearly-flat, 
"Mexican-hat-shaped" mirror configuration, K. S. Thorne proposed a nearly- concentric mirror con- 
figuration capable of producing the same mesa beam profile on the mirror surfaces. Subsequently, 
Bondarescu and Thorne introduced a generalized construction that leads to a one-parameter fam- 
ily of "hyperboloidal" beams which allows continuous spanning from the nearly-flat to the nearly- 
concentric mesa beam configurations. This paper is concerned with a study of the analytic structure 
of the above family of hyperboloidal beams. Capitalizing on certain results from the applied optics 
literature on flat-top beams, a physically-insightful and computationally-effective representation is 
derived in terms of rapidly-converging Gauss-Laguerre expansions. Moreover, the functional relation 
between two generic hyperboloidal beams is investigated. This leads to a generalization (involving 
fractional Fourier transform operators of complex order) of some recently discovered duality rela- 
tions between the nearly-flat and nearly-concentric mesa configurations. Possible implications and 
perspectives for the advanced LIGO optical cavity design are discussed. 

PACS numbers: 04.80.Cc, 07.60.Ly, 41.85.Ew, 42.55.-f 



I. INTRODUCTION 

The current baseline design for the Laser Interferom- 
eter Gravitational-wave Observatory (LIGO) Q, as well 
as that for its advanced version , is based on the use of 
Fabry-Perot optical cavities composed of spherical mir- 
rors, which support standard Gaussian beams (GBs) 1^. 
During the past few years, there has been a growing 
interest toward the use of non-spherical mirrors as a 
possible aid for reducing the thermal noise of the mir- 
rors 4J. In particular, it was proposed by D'Ambrosio, 
O'Shaughnessy and Thorne to replace the GB profile 
with a flat-top (commonly referred to as "mesa") profile, 
for better averaging the thermally-induced mirror surface 
fluctuations. They showed that such mesa beams could 
be synthesized via coherent superposition of minimum- 
spreading GBs with parallel optical axes, and could be 
supported by nearly-flat, "Mexican-hat-shaped" mirrors 
|38|. A thorough investigation of the theoretical impli- 
cations and implementation-related issues |E S IE Si EOI 
indicated a potential reduction by a factor three in the 
thermoelastic noise power and a factor two in the coating 
Brownian thermal noise power, without substantial fab- 
rication impediments. A prototype optical cavity is cur- 
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rentl y b eing developed, and experimental tests are under 
way 

In gravitational-wave interferometers, a serious con- 
cern is posed by the light-prcssure-induced tilt-instability 
of the cavity mirrors. In this connection, the inherent 
tilt-instability of the current (nearly-flat, spherical mir- 
ror) LIGO baseline design, was flrst pointed out by Si- 
dles and Sigg |3- Subsequently, Savov and Vyatchanin 
[l3| found similar effects for the nearly-flat mesa (FM) 
designs. Based on these observations, and on the re- 
sults by Sidles and Sigg concerning the compar- 
ison between nearly-flat and nearly-concentric spher- 
ical mirrors, Thorne proposed an alternative nearly- 
concentric Mexican-hat-shaped mirror configuration, ca- 
pable of supporting mesa beams with intensity distribu- 
tion at the mirror identical with that of the FM configura- 
tion, but featuring a much weaker j39| tilt-instability |T3| . 
These nearly-concentric mesa (CM) beams are synthe- 
sized by coherent superposition of minimum-spreading 
GBs with non-parallel optical axes sharing a common 
point. Quite remarkably, the FM and CM configura- 
tions were found to be connected through a duality re- 
lation, first discovered numerically by Savov and Vy- 
atchanin p^ . and subsequently proved analytically by 
Agresti et al. [T^ . which allows a one-to-one mapping 
between all the corresponding eigenmodes. The geomet- 
rical construction underlying FM and CM beams was fur- 
ther generalized by Bondarescu and Thorne , in terms 
of a family of "hyperboloidal" beams, parameterized by 
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a "twist-angle" a G [0, tt] which allows continuous span- 
ning from the FM (a = 0) to the CM (a — tt) configura- 
tions, passing throu gh t he standard GB (a = 7r/2) case. 
It was suggested in that the optimal configuration, 
in terms of both thermal-noise and tilt-instability reduc- 
tion, should be found in a neighborhood of a = tt (CM 
configuration). This renders the family of Bondarescu- 
Thorne (BT) hyperboloidal beams of potential interest 
in the design of advanced LIGO. 



This paper elaborates on the analytic structure of the 
BT hyperboloidal beams. Our investigation capitalizes 
on and generalizes a number of results from the applied 
optics literature concerning flat-top beams, which have 
most likely not come to the attention of the gravita- 
tional wave community. Indeed, during the past decade, 
flat-top beams have drawn a considerable attention from 
the applied optics community, and several models have 
been proposed and investigated. Prominent among them 
are the celebrated "supergaussian" beams 3,, J£j, the 
"flattened Gaussian" beams introduced by Gori and co- 
workers 0, 0, the "flattened" beams introduced 
by Sheppard and Saghafi 20J, the "flat-topped multi- 
Gaussian" beams introduced by Tov ar [21[ . and the "flat- 
topped" beams introduced by Li |22L |23|| [40l |. With 
the exception of the analytically-intractable supergaus- 
sian, all other models admit analytic parameterizations in 
terms of Gaussian , Gaussian-Laguerre (GL) 

[UllliOl, or "elegant" GL jlil beam expansions. In 
this paper, we first show that the FM and CM beams 
belong to the class of flattened beams introduced in [23| , 
and can therefore be represented in terms of the rapidly- 
converging GL beam expansions derived therein. Based 
on this observation, we then generalize the approach in 
[20| to accommodate the more general family of BT hy- 
perboloidal beams j^l- This leads to a generalization 
(at least for the dominant eigenmode) of the duality rela- 
tions discovered in |l4| , which involves fractional Fourier 
transforms of complex order. The above results, here 
discussed for the simplest case of the dominant eigen- 
mode, set the stage for the development of new problem- 
matched computational tools for the modal analysis of 
Fabry-Perot optical cavities supporting general BT hy- 
perboloidal beams. 



The remainder of the paper is laid out as follows. Sec- 
tion ^] introduces the problem geometry, and provides a 
compact review of background results from 15] on BT 
hyperboloidal beams and supporting mirrors. Section UTTI 
contains the analytic derivations concerning the GL ex- 
pansions and the generalized duality relations, as well as 
representative numerical results for validation and cali- 
bration. Section llVl contains preliminary conclusions and 
recommendations. 



nearly-spheroidal mirrors 
waist plane 




FIG. 1: Problem schematic: A perfectly symmetric Fabry- 
Perot optical cavity composed of two nearly-spheroidal mir- 
rors separated by a distance L along the z-axis. The trans- 
verse coordinates at the waist {z — 0) and mirror {z = L/2) 
planes are denoted by ro and r, respectively. For the a- 
parameterized family of hyperboloidal beams of interest, the 
mirror shape is obtained by adding the perturbation —ha in 
to the fiducial spheroid in 



II. BACKGROUND: BT HYPERBOLOIDAL 
BEAMS AND SUPPORTING MIRRORS 

In this Section, we briefly review the procedure pro- 
posed by Bondarescu and Thorne [Tsl l for constructing a 
family of "hyperboloidal" beams, which contains as spe- 
cial limiting cases the FM and CM beams. Referring to 
the problem geometry illustrated in Fig. ^ we consider a 
perfectly symmetric Fabry-Perot optical cavity composed 
of two nearly-spheroidal mirrors separated by a distance 
L along the z-axis of a Cartesian {x, y, z) (and associated 
cylindrical (r, 6*, z)) coordinate system. The transverse 
coordinates at the waist {z = 0) and mirror (z = L/2) 
planes are denoted by Tq = xqx -I- yoY = ''o cos O^vi -\- 
rosin^oy and r = xx + yy — rcos^x -I- rsinOy, re- 
spectively. Here and henceforth, x, y and z denote the 
standard Cartesian unit vectors. Throughout the paper, 
an implicit time-harmonic exp(— io^i) dependence is as- 
sumed for all field quantities. 

The BT construction [l5], which generalizes the origi- 
nal idea in ^ , is based on the superposition of minimum- 
spreading GBs launched from a circular equivalent aper- 
ture of radius Rq at the waist plane (z — 0), with optical 
axes pointing along the unit vector 

Uq = — [cos 6*0 — cos(6'o — Oi)] x 
L 

-\- J- [sin 6*0 - sin(6'o ~a)]y + z. (1) 

As shown in Fig. |2 these optical axes are the genera- 
tors of hyperboloids. The "twist-angle" a in param- 
eterizes this family of "hyperboloidal" beams, allowing 
continuous spanning from the FM configuration (a = 0, 
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(a) 




In lO, A is an a-independent complex constant, and 
Wq is the GB spot size at the waist. According to the 
minimum-spreading criterion, this spot size is chosen as 



(4) 



"^0 = \ -r- 



where fco = 27r/Ao denotes the free-space wavenumber 
(Aq denoting the free-space wavelength), so that the mir- 
ror plane is located exactly at the Rayleigh distance 0, 



(b) 



(c) 




FIG. 2; Geometrical construction of the BT hyperboloidal 
beams in (|3J: Minimum-spreading GBs are launched from a 
circular equivalent aperture of radius Ro at the waist plane 
{z ~ 0), with optical axes pointing along the unit vector Ua 
in 10. (a) FM beam (a = 0): Optical axes are the generators 
of a cylinder, (b) CM beams (a = n): Optical axes are 
the generators of a cone, (c) Generic hyperboloidal beams 
[0 < a < n): Optical axes are the generators of a hyperboloid. 



cylindrical degenerate, cf. Fig. Ufa)) to the CM configu- 
ration {a — TT, conical degenerate, cf. Fig. Elb)). For this 
family of beams, the wavefronts at the mirror location are 
roughly approximated by the "fiducial" spheroids 
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sin^(Q!/2) 



r2 sin2(a/2) 



, r«L/2, 



(2) 



which degenerate into planar and spherical surfaces in the 
FM la = 0) and CM {a = tt) case, respectively. Follow- 
ing [13, an integral expression (valid under the paraxial 
approximation) for the (unnormalized) beam field distri- 
bution on these fiducial surfaces can be written as 



Ua{r,Sa) = A / dro / (i6'oro exp 

JQ Jo 



rro . . 

I — 2" sm Uq sma 



(r^ +rl- 2rrocos6'o) . , 

-^^ —2 (1 — « cos a) 



(3) 



ZR 



kowl 



L 
2' 



(5) 



Note that the expression in ||2Jl is valid only on the fiducial 
surface z = Sa{r). For a — 7r/2, the double integral in 
Q can be computed in closed form, yielding a simple 
Gaussian lisll. 



(6) 



= Aocxp ( 





= C/a, 


U-^ — a 




A 


A* 



with Aq denoting a complex constant. For other values of 
a, the radial integral in (PJ can still be computed analyt- 
ically, whereas the angular integral has to be evaluated 
numerically. It is readily verified from (|3Jl that the fol- 
lowing symmetry relations hold: 



(7a) 
(7b) 



where * denotes complex conjugation; this sets the mini- 
mal meaningful range for the twist- angle a to [0, tt]. The 
relation (|7b|) can be interpreted in the broader duality 
framework detailed in (see Section Fill A 111 . 

From the theory of graded-phase mirrors ^4, 25], it is 
well-known that, in order for an optical cavity to sup- 
port a stable beam with a given profile as the fundamen- 
tal eigenmode, its mirror profile has to match the beam 
wavefront. For the BT hyperboloidal beams in this 
can be achieved by applying a correction 

arg[f7„(r,5„)] - arg[[/„(0,S'a)] 
ha(r) = ; 8 



to the fiducial spheroidal shape Sa in (0), so that 
aTg[Ua{r, Sa — ha)] = constant. 



(9) 



For the FM (a = 0) and CM (a = tt) cases, the correction 
in JHl) reduces to the Mexican- hat-shaped profile in 0] 
(see also Fig. Ela) below). Moreover, from Ij7bp . the 
remarkable result 



K~a{r) = -ha{r) 



(10) 



follows, which can also be interpreted within the above- 
mentioned duality framework 



erpreted 



4 



III. ANALYTIC STRUCTURE OF BT 
HYPERBOLOIDAL BEAMS 

Capitalizing on the background results summarized in 
Section |nl in this Section, we develop the analytic GL 
representation of general BT hyperboloidal beams (valid 
at any point in space, within the limit of the paraxial ap- 
proximation) as well as of the supporting mirror profiles. 
Moreover, we generalize the symmetry/duality relations 
in lO and (|10() to the most general case involving two 
arbitrary values of the twist-angle parameter. For such 
generalization, we provide functional and optical inter- 
pretations, based on fractional Fourier operators of com- 
plex order. 



A. Field Distribution at Waist 

We begin by considering the field distribution at the 
waist plane (z = 0) for the FM (a = 0) beam 0,[l3|, 



Uo{r,0) 



drQ exp 



■ro<-Ro 



i-ol 



2 f"" /2rro 
— / droro/o — I'^^P 



"'0 



(r^ + rjp 



, (11) 



where Iq (f ) denotes a zeroth-order modified Bessel func- 
tion of the first kind [H Sec. 9.6]. 



1. Duality Relations 

In p^ . within a broader framework of duality rela- 
tions, the CM {a — tt) beam field distribution at the 
waist plane was shown to be related to the FM one in 
via a Fourier transform operator 



n 



[/,(r,0) ^ ;7o(r,0) 



(12) 



which, in view of the assumed cylindrical rotational sym- 
metry, takes the form of the Hankel transform (HT) 

[F{r)] = \ [ droroF{ro)Jo (^) • (13) 
Wq Jo \ J 

In (|13|l and henceforth, J„ (£) denotes an nth-order Bessel 
function of the first kind ^9, Sec. 9.1]. Straightforward 
application of the HT l(T^ to ((TT|l yields, via the convo- 
lution theorem 



U^{r,0) = [Uo{r,0)] 
rRo ^ V f^o 



exp 



-2 • (14) 



2. GL Representation 

The CM field distribution in H14|l is recognized to coin- 
cide with the one used in ,20, Sec. 2] to generate beams 



with a "flattened" far-field profile, for which a GL beam 
expansion was subsequently derived. Following ^20] . the 
field distribution in H14|l can be expanded as 



C/.(r,0)= ^A^Va 



' V2r 

Wq 



(15) 



In p5|l . tpmiO ^''6 orthonormal GL basis functions, 

V^„(^) = x/2exp f-^\ L,n{e), (16a) 



(16b) 



where Ln{Q denotes an nth-order Laguerre polynomial 
[2^ Chap. 22], and Spq denotes the Kronecker symbol. 
The expansion coefficients Am-' in (|15|l are given by 



%/2wg 
Rl 



Rl 



^ "^ + l':dkh (17) 



2wl 



where P(n, ^) denotes an incomplete Gamma function 
|26L Eq. (6.5.13)]. The behavior of the expansion co- 
efficients in (|17|l . as a function of the summation in- 
dex m, is shown in Fig. for three representative 
values of the ratio wq/Rq: The Am"* are almost con- 
stant for m < _Rn/(2wo), and fall off quite abruptly for 
TO > Rq/{2wq) 41j. For the parametric range of interest 
for the LIGO design {wq/Rq = 0.25), this results in a 
rapidly-converging (m < 20) expansion H15I) . 

The corresponding GL expansion for the FM field dis- 
tribution at the waist plane can be obtained by exploiting 
the HT relation in (|12|l . Similar arguments were likewise 
invoked in [23, Sec. 4], in the form of near-to- far-field 
transformations, to generate beams that were flattened 
at the waist plane ("inverted flattened" beams, in the 
notation of ^^)- In this framework, one first observes 
that the GL basis functions in H16a|l are eigenfunctions 
of the HT operator in ifBjl (271] . 



H 



Wo 



4>r, 



V2r 

Wo 



(-l)"V'r, 



V2r 

Wq 



(18) 



Equation H18|) can be derived as a special case (corre- 
sponding, e.g., to letting a = 1, 6 = 2, ^ = 2r/wo, 
( = ro/wo) of the identity [H p. 43] 



Cexp(-aC')i„(feC')Jo(CC)dC 



be 



Aa{b - a) 



_ (a - 6)™ 
^ 2a"'+i 

, Re(a) > 0. (19) 



Application of the HT to then yields, via l(TF|l . 
the corresponding GL expansion for the FM field distri- 
bution in Hll|l . with the following mapping between the 
expansion coefficients; 



= (-l)"yl('"). 



(20) 



5 



0.01 



lE-3 




lE-4 



lE-5 



FIG. 3: GL expansion coefficients in 11711 vs. summation 
index m, for different values of wo/Ro- Continuous curve: 
Wo/Ro = 0.5; Dashed curve: wq/Rq = 0.25; Dotted curve: 
wo/Ro = 0.1. 



Alternatively, as shown in Appendix ^ one could derive 

directly the FM expansion coefficients Am'' , and then ex- 
ploit the mapping in (|20|l to obtain those pertaining to 
the CM configuration. The above derivations clarify the 
relationship between the CM and FM beams in 
and the GL expansions for the "flattened" and "inverted 
flattened" beams in |23| , respectively. A natural question 
then arises, as to whether the expansion coefRcient map- 
ping in H20() can be generalized to arbitrary values of the 
twist-angle a, thereby allowing a GL representation for 
general BT hyperboloidal beams. Recalling that, from 



0, TO > 0, 



one notes that the mapping 

AL"' = (- cos ar^W 



(21) 



(22) 



accounts correctly for the three notable cases a = 
(FM), a = 7r/2 (GB) and a = tt (CM). One is accord- 
ingly led to speculate whether the GL expansion 



c/„(r,o) = 



m— 



' V2r' 

Wo 



(23) 



may hold for arbitrary values of the twist angle a. This 
turns out to be indeed the checked by numerical 

comparison against the BT reference solution in (see 
Section ImPl below). The analytic GL expansion in l|23|l 
is obtained here for the first time, to the best of our 
knowledge, and sets the stage for a generalization, to 
arbitrary values of the twist-angle, of the duality relation 
in H12|l . whose possible interpretations and implications 
are discussed below. 



3. Generalized Duality Relations: Functional and Optical 
Interpretations 

We begin by considering a class of cr-parameterized 
modified HT operators defined as 



(1 + cr) Jo 



rodroF{ro)Jo 



xexp 



{r^ + rl){l-<y) 



wl{l + cj) 



a > 



irro^/a 
-1. (24a) 



For (T < — 1, the integral in (|24a|l diverges for the beams 
of interest here (decaying as 0[exp(— t^/wq)] in the waist 
plane), and the following definition should be used: 

ni^oHFir)]^ni-^^ {H«[i^(r)]} 

= <){hL7'[^M]}, <^<-l- (24b) 

The operator in H24|) generalizes the ordinary HT in 
(|13|l : it is readily verified that it reduces to the ordinary 
HT for cr — 1, and to the identity operator for a — —1. 
From (|19() . it then follows the generalization of the 
eigenproblem in (|18|) . 



^/2r 

Wo 



= (-^)™V^„ 



V2r 

Wo 



(25) 



Application of the generalized HT l|24|l to the GL ex- 
pansion in l|23|l reveals, via (|25() . the functional relation 
between the field distributions at the waist plane per- 
taining to two BT hyperboloidal beams characterized by 
generic values, ai and a2, of the twist-angle. 



[/„,(r,0) ^ ;7„,(r,0). 



cosa2 
cos ai 



(26) 



The generalized HT in H26() extends the duality relation 
in ifT^ to the most general case, and admits a suggestive 
analytic interpretation in terms of (the cylindrical version 
of) a fractional Fourier operator of complex order |2^ 
Isfi Isil 1^3 ■ From the physical point of view, complex- 
order Fourier transform operators can be interpreted in 
terms of propagation through a paraxial optical system 



described "29' 
i43j. For ^ 



'30|, |31|, |33| by a complex ABCD matrix 
the ABCD matrix can be shown to be 



,.(1 



(1 + ^) 
kowl^ 
/7r7\ 



fcowg(l + cr) ' 
(1-a) 



cos 
2 

sm 

. L 



2) 



2V^ 

L /7r7 
— sm — 
2 V 2 

cos(- 



(27a) 



(27b) 



In Ij27b|l . the parameter 7 denotes the complex order of 
the transform, defined as l23| 



.2, 

-I- log 



.log(cr) 



(28) 
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Mr) 



Im{r) 



a, 
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= 1 
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> )•( 
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7r/2 









FIG. 4: Admissible values of the complex order 7 1281 1 of the 
generalized HT operator in (I24II relating the field distributions 
of two BT hyperboloidal beams with generic twist-angles, qi 
and Q2 (cf. l^bU 'l. The real part of 7 can either be 1 or 
0, whereas the imaginary part is generally nonzero, except 
along the two diagonals a\ = 02 (o" = —1, 7 = 0) and ai = 
71 — 0:2 ((7 = 1,7 = 1), where the generalized HT operator 
in 12411 reduces to the identity and ordinary HT operator, 
respectively. 



The admissible values of the complex order in (|28|l are 
illustrated schematically in Fig. 21 over the meaning- 
ful (ai, a2)-range. It is observed that the real part can 
either be 1 or 0, whereas the imaginary part is gener- 
ally nonzero, except along the two diagonals ai = a2 
and ai — tt — a2 where the generalized HT operator in 
(|24|l reduces to the identity and ordinary HT operator, 
respectively. 

We are currently exploring the possibility of extending 
the above results, limited to the dominant cigcnmodc, to 
higher-order (isotropic and nonisotropic) modes. If suc- 
cessful, such an extension would provide valuable phys- 
ical insights and computational savings, unraveling the 
effects of the twist-angle parameter in the eigenspectrum 
(see also the discussion in Section llV|) . 



B. Paraxial Field Distribution: GL Beam 
Representation 

From the GL expansion at the waist plane in (|23l) . a 
paraxially-approximated expression for the field distribu- 
tion of a generic BT hyperboloidal beam at any point in 
space can readily be obtained as 



(29) 



rn— 



where {r, z) denote the standard GL beam propaga- 
tors 



w{z) 



Ipr, 



V2r 



w{z) 



exp 



2R{z) 



X exp{i[/co^- (2m-f l)$(z)]}, (30) 



which match the GL basis functions at the waist plane, 
^^{r,0) = i:^{V2r/wo)- In w{z), R{z) and $(z) 

denote the standard GB spot size, wavefront radius of 
curvature, and Gouy phase, respectively 



w{z) = wqx 1 + 



- , i?(z) = z + ^, 
zrJ z 



$(z) = arctan ( — 



(31) 



with the Rayleigh distance zn defined in 0. The GL 
beam expansion in H29() represents, together with the gen- 
eralized duality relation in (|26l) . the main original result 
in this paper. We stress that, unlike the expression in ((^J, 
the representation in ()29|l is valid at any point in space, 
within the limits of the paraxial approximation. When 
evaluated on the fiducial surface z = Sa (r) , it yields 



(1 - i) 
Ua{r,Sa) ~ — - — exp 



, L r^cosa 



00 / 



m— 



(32) 



where the approximation 5*^ ~ i/2 has been used in 
(|31|l . It is easily verified that (|32|l satisfies the phase- 
conjugation symmetry relation in l|7b|l . Moreover, by 
comparing (|32|l to (|23|l , and recalling (|26|l , one can easily 
derive a generalized duality relation for the field distribu- 
tions on the fiducial surfaces, in terms of the generalized 
(complex-order) HT in (|^ . 



n 



exp 



kor cos a2 
2L 



exp 



fcor'^cos ai 



2L 



cos a2 
cos ai 



(33) 



The relation in (|33|l (which is similar to that in l|26|l . 
apart for the phase factors and a scaling by a factor \/2 
in the GB spot size) generalizes completely the symmetry 
relations in {Tj). 



Mirror Profile 



The correction to be applied to the fiducial spheroidal 
mirror shape Sa in (0) is finally obtained by substituting 
into 



ha{r)K— arg 



expi I 



fegr^cos a 
~2L 



00 ^ V 



(34) 
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It is readily verified that, in view of H18|l . the correction 
profile in (|34|l satisfies the duality relation in (|10|) . Trun- 
cation of the infinite series in H34(l is not an issue, as 
further discussed below. 



D. Representative Results 

In order to validate and calibrate the proposed GL rep- 
resentations, we now move on to illustrating some rep- 
resentative numerical results. In all examples below, all 
the relevant parameters were chosen as in 0, Il5j . More 
specifically: L — 4km (length of the optical cavity, cf. 
Fig. Ao = 1064nm (wavelength of the laser beam), 
Wo = ^JL\q/{2'k) — 2.603cm (GB spot size at waist), 
i?o — 4u'o = 10.4cm (radius of the equivalent aperture 
distribution at the waist plane). For the truncation of 
the GL series involved, a simple criterion was utilized, 
requiring that the magnitude of the last retained M-th 
term is less than 0.1% of that of the leading term. 



A 



(a) 
M 



A. 



(a) 



< 10" 



(35) 



For the cases a = 0,7r (see the dashed curve in Fig. |2J), 
this yields M — 18. In view of the coefficient mapping 
in H20|l . the convergence becomes faster as a approaches 
the critical value of ■n/2 (pure GB, for which one obtains 
only one nonzero coefficient). 

As a reference solution, we considered the BT inte- 
gral representation in (O, where the radial integral was 
computed analytically, and the angular integration was 
performed numerically utilizing the adaptive quadrature 
routines of Mathematica® 33]. 

Some representative results for the field distribution 
are shown in Fig. [3 Specifically, Fig. EJa) shows the GL- 
computed (via lj23)) intensity distribution on the fiducial 
surface, for various values of the twist-angle a, illustrat- 
ing the gradual transition from Gaussian (a = 7r/2) to 
mesa (a = 0,7r) profile. To quantify the agreement with 
the reference solution, Fig. El^b) shows the relative error 



(36) 



where the superfix ^^^^^ denotes the BT representation 
in (where the complex constant A is determined by 
enforcing the matching with the GL expansion at r = 0). 
The error (|36|l never exceeds 0.1%, over the region of sig- 
nificant field intensity (and drops below numerical preci- 
sion for the a = tt/2 pure GB case). This is consistent 
with the truncation criterion in (|35|l . which can therefore 
be used to control the accuracy. 

The results pertaining to the mirror profiles are shown 
Fig. El Specifically, Fig. EJa) shows the corrections 
ha computed via (|34|) . illustrating the gradual transi- 
tion from the spherical (a — tt/2) to the Mexican- hat 
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FIG. 5: BT hyperboloidal beam field distribution evaluated 
on the fiducial surface z — Sa{r), for diflterent values of the 
twist-angle parameter a. Optical ca vity parameters: L = 
4km, Ao = 1064nm, wq = VLAo/(27r) = 2.603cm , and Ro = 
Awo = 10.4cm. (a) Intensity distribution computed via 11321 . 
using the truncation criterion in 13511 . (b) Relative error in 
13511 . Continuous curve: a — 0,-k{M = 18); Dashed curve 
a = O.Itt, 0.97r(M = 17); Dotted curve: a = 0.2tv, 0.8tt{M = 
14); Dotted-dashed curve: a = 0.5tv{M = 0). 



(a = 0, 7r) mirror profile. Figure|SJb) shows the absolute 
error 



(37) 



with respect to the reference solution. The error 13 7|) 
never exceeds 10~^Ao over the significantly illuminated 
portion of the mirror. For the LIGO design (Ao = 
1064nm), this corresponds to errors ^O.lnm, well within 
the typical fabrication tolerances. 



IV. CONCLUSIONS AND 
RECOMMENDATIONS 

In this paper, the analytic structure of a family of 
hyperboloidal beams, introduced by Bondarescu and 
Thorne [l5l | as a generalization of the mesa beams sup- 
ported by Mexican-hat-shaped mirrors, has been inves- 
tigated. Rapidly converging expansions in terms of GL 
beams have been first introduced for the "extremal" cases 
of FM and CM beams, capitalizing on results from [20j . 
The representation has been then extended to the more 
general BT hyperboloidal beams in ^l^], leading to a 
complete generalization (for the dominant eigenmode) of 
the duality relations introduced in |l4j . based on frac- 
tional Fourier transforms of complex order. The above re- 
sults, numerically validated and calibrated against a ref- 
erence solution independently-generated from 'l5|, pro- 
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FIG. 6: Nearly-spheroidal mirror profiles supporting the BT 
hyperboloidal beams in Fig. |3 (a) Correction ha to the fidu- 
cial spheroid Sa in computed via I34II . using the trunca- 
tion criterion in l|35|l . (b) Absolute error in 1)37^ . Continuous 
curve: a = n{M = 18); Dashed curve: a = 0.97r(M = 17); 
Dotted curve: a — 0.8tt{AI = 14); Dotted-dashed curve: 
a = 0.57r(A/ — 0). The correction profiles pertaining to 
a = 0, O.Itt, 0.2tt (not shown) differ merely by sign from those 
pertaining to a = tt, O.Qtt, O.Stt, respectively (cf. I)l()|l ). 



vide a physically-insightful and computationally-efFective 
parameterization of the beam and mirror profiles. It is 
hoped that they may help addressing the optimization 
of the advanced LIGO optical cavities in a broader per- 
spective. In this framework, current and future research 
directions include: 

i) Thorough parametric analysis of the family of BT 
hyperboloidal beams, as well as of other classes 
of flat-top beams Ulllllllllllllilllllll, 

aimed at finding optimal design criteria in terms of 
thermal-noise and tilt-instability reduction. 

ii) Development of semi-analytic, problem-matched 
techniques for the computation of higher-order 
eigenmodes in nearly-spheroidal-mirror optical cav- 
ities supporting general BT hyperboloidal beams. 
These techniques should take advantage from global 
GL expansions, as compared to local discretization 
schemes presently in use H, 01 • 

Hi) Full extension of the duality relations in to 
the family of BT hyperboloidal beams and support- 
ing mirrors. Such an extension, demonstrated here 
for the dominant eigenmode, should be based on 



the "complexification" of the order of the involved 
Fourier transform operators. Finding such a one- 
to-one mapping between eigenmodes with arbitrary 
values of the twist-angle parameter would provide 
important physical insight and computational ad- 
vantages. 
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APPENDIX A: ALTERNATIVE DERIVATION OF 
FM GL EXPANSION COEFFICIENTS 

It is instructive to illustrate the direct derivation of the 
GL expansion coefficients pertaining to the FM configu- 
ration, starting from the integral representation in 
In this framework, one first expands the modified Bessel 
function Iq in series of Laguerre polynomials, 



(Al) 



where ^ — \/2?'/wo, and 



!^ L„(e2)exp(-C')/o|^ 



(-1)" / 



, 2m 



ro 



l\ 22'^ 



Wo 



Wo 
„2 



Cd^ (A2a) 



exp(^). (A2b) 



Equation ljA2ap follows from the Laguerre polynomials 
orthogonality condition 26, Eq. 22.2.13], whereas (jA2b|) 
can be derived from .35;, Eq. (2.19.1.17)]. Substituting 
(|A1|) and (|A2bp into (|ll|l then yields, after some manip- 
ulation, the GL representation 

°° /n 2 



m=0 



X / exp(-C)rdC 
/o 



(A3) 



whose expansion coefficients, recalling the integral ex- 
pression of the incomplete Gamma function |2a . Eq. 
(6.5.1)], coincide with those in (|20|) . 
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